Phase of bi-particle localized states for the Cooper problem in two-dimensional 

disordered systems 
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The Cooper problem is studied numerically for the Anderson model with disorder in two-dimensions. 
It is shown that the attractive Hubbard interaction creates a phase of bi-particle localized states in 
the regime where non-interacting states are delocalized. This phase cannot be obtained in the mean- 
field approximation and the pair coupling energy is strongly enhanced in this regime. The effects of 
magnetic field are studied and it is shown that under certain conditions they lead to derealization. 

PACS numbers: 74.20-z, 74.25-q, 74.40+k 

I. INTRODUCTION 

Recently a great deal of attention has been attracted 
to investigation of superconductor-insulator transition 
(SIT) in systems with disorder. Various approaches are 
used to study this problem including analytical theo- 
retical methods Jl],^], intensive numerical simulations of 
many-body quantum systems with quantum Monte Carlo 
methods ]3|,|| , as well as mean field numerical simulations 
These theoretical studies are stimulated by chal- 
lenging experiments on SIT in disordered films and 
high-T c superconductors |Ic|Jli|] . The results obtained 
in Jl0| show an interesting correlation between the opti- 
mal doping and the Anderson transition in the normal 
phase obtained by application of a strong pulsed mag- 
netic field. Even if the experiments fl(]|ll|] arc done 
with three-dimensional crystals the coupling between 
two-dimensional planes is relatively weak and the two- 
dimensional effects should play an important role. Due 
to that it is relevant to study the SIT in two-dimensional 
disordered systems. In the case of weak disorder the 
Anderson theorem Jl^,|l3| guarantees that the supercon- 
ductivity is not affected by disorder. However it is not 
obvious if the theorem is still valid in the presence of rel- 
atively strong disorder. It is quite possible that in this 
regime the interplay of disorder and interaction can lead 
to appearance of new physical effects. The theoretical in- 
vestigation of this regime is however rather difficult. The 
existent analytical methods are not well adapted to the 
regime of strong interaction and disorder. At the same 
time the numerical studies also meet serious difficulties. 
Indeed the direct diagonalization methods are restricted 
to relatively small system size since the Hilbert space 
grows exponentially with the number of particles Jl4],|l5| . 
The quantum Monte Carlo methods are not so sensitive 
to a huge size of the Hilbert space but still they are re- 
stricted to systems of quite moderate size (for example 
lattices of 8x8 sites in 

In the view of above numerical difficulties it is natu- 
ral to develop the approach introduced by Cooper |lfj| ] 
and to study the problem of two particles with attrac- 



tive interaction near the frozen Fermi sea in the pres- 
ence of disorder. Even if the original Cooper problem 
of two particles without disorder does not reproduce ex- 
actly the BCS theory of many-body problem it never- 
theless captures the essential physical properties of the 
system and gives the appearance of coupled states with 
a qualitatively correct coupling energy and correlation 
length. Without disorder the Cooper problem can be 
solved exactly. However in presence of disorder the situ- 
ation becomes more complicated even for only two inter- 
acting particles (TIP). Indeed, for relatively strong dis- 
order even the matrix elements of interaction between 
non-interacting eigenstates can not be obtained analyt- 
ically and due to that the problem should be studied 
numerically. The first numerical studies of the Cooper 
problem in the presence of disorder were done in |T^ ] 
for two particles with attractive Hubbard interaction in 
the three-dimensional Anderson model. These studies 
showed that the interaction can lead to localization of 
pairs in the non-interacting metallic phase. This result 
is qualitatively different from the mean field solution of 
the Cooper problem in the presence of disorder (Cooper 
ansatz) which gives delocalized pairs for the same pa- 
rameters. This shows that the non-diagonal interaction 
induced matrix elements play an important role and lead 
to new physical effects which are not captured by the 
mean field approximation. 

In this paper we study the Cooper problem on a two- 
dimensional lattice with disorder described by the An- 
derson model. Our numerical studies show that near 
the Fermi level the attractive Hubbard interaction be- 
tween two particles creates localized pairs in the regime 
where non-interacting eigenstates are well delocalized 
(extended). The coupling energy of these pairs is much 
larger than the coupling energy given by the mean field 
solution (Cooper ansatz). Therefore energetically it is 
more favorable to have an insulator with localized pairs 
instead of usual weakly coupled delocalized Cooper pairs. 
This result indicates the appearance of a new phase of bi- 
particle localized states (BLS phase) which appears in the 
regime when non-interacting states are extended (metal- 
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lie). It is in a qualitative agreement with the quantum 
Monte Carlo studies obtained recently in Toulouse [fTs) . 
This BLS phase is qualitatively different from BCS so- 
lution which corresponds to weakly coupled delocalized 
pairs. The paper is organized as follows. The proper- 
ties of BLS phase without magnetic field are discussed in 
Section 2. The effects of perpendicular magnetic field on 
the ground state properties in the presence of interaction 
and disorder are analyzed in Section 3. The discussion 
of the results is presented in the last Section. 

II. GROUND STATE PROPERTIES WITHOUT 
MAGNETIC FIELD 

To study the Cooper problem of two interacting 
particles in the presence of disorder we use the two- 
dimensional Anderson model. In this model the one- 
particle eigenstates are determined by the Hamiltonian 

H l = J2E Il \n){n\ + V £ |n)(m| (1) 

n <n,m) 

where n and m are index vectors on the two-dimensional 
square lattice with periodic boundary conditions, V is 
the nearest neighbor hopping term and the random on- 
site energies E a are homogeneously distributed in the en- 
ergy interval [— -y-, , where W is the disorder strength. 
This one-particle model has been extensively studied by 
different authors, see for example [jl9). For the two- 
particle problem on this two-dimensional lattice we con- 
sider on-site attractive Hubbard interaction between par- 
ticles of strength U < 0. We consider the particles in 
the singlet state with zero total spin so that the spatial 
wavefunction is symmetric with respect to particle per- 
mutation (interaction is absent in the triplet state). 

To investigate the effects of interaction between par- 
ticles near the Fermi level we generalize the Cooper 
approach for the case with disorder. To do that we 
rewrite the TIP Hamiltonian in the basis of one-particle 
eigenstates of the Hamiltonian ([[]). In this basis the 
Schrodinger equation for TIP reads 

{E nll + E m2 )x mi ,m 2 + U Qm 1 ,m 2 ,m' 1 ,m' 2 Xm' 1 ,m' 2 

m 1 ,m 2 

= EXmi,m 2 - ( 2 ) 

Here E m are the one-particle eigenenergies correspond- 
ing to the one-particle eigenstates \4> m ) and x mi , m2 
are the components of the TIP eigenstate in the non- 
interacting eigenbasis \4> mi , 4>m 2 )- The matrix elements 
UQ mi m2 m ' m §i yes the interaction induced transi- 
tions between non-interactive eigenstates \4> mi , 4>m 2 ) an< i 
\d> ' . d> '). These matrix elements are obtained by 
rewriting the Hubbard interaction in the non-interactive 




FIG. 1. Ground state probability distributions of two in- 
teracting particles for the Cooper problem with disorder on 
a lattice of linear size L = 40. The cases (a,b) show one par- 
ticle probability distribution /(n) in absence of interaction 
(U = 0V) for the disorder strength W = 2V (a) and W = 5V 
(b). All other cases are obtained for the Hubbard interaction 
U = — 2V. The cases (c,d) show the one-particle probability 
/(n) (c) and the interparticle distance probability /d(r) (d) 
for W = 2V . The same probabilities are shown in cases (e,f) 
for W = W (/(n) for (e) and / d (r) for (f)). The cases (g,h) 
present the probabilities for the same W = 5V as in cases 
(e,f) however here the ground state is obtained in the mean 
field approximation of the Cooper ansatz. All data are given 
for the same realization of disorder. 

eigenbasis of model (nh. In the analogy with the origi- 
nal Cooper problem |16| | the summation in (j^) is done 
over the states above the Fermi level with eigenenergies 
E m > > Ep with m 12 > 0. The Fermi energy Ep w is 

determined by a fixed filling factor v = 1/2. To keep the 
similarity with the Cooper problem we restrict the sum- 
mation on m 1 2 by the condition 1 < m 1 + m 2 < M. In 
this way the cut-off with M unperturbed orbitals intro- 
duces an effective phonon frequency ujd oc Mj L 2 = 1/a 
where L is the linear system size. When varying L we 
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keep a fixed so that the phonon frequency is independent 
of system size. All the data in this work are obtained 
with a = 15 but we also checked that the results are not 
sensitive to the change of a. We note that a similar TIP 
model was considered for the problem of two repulsive 
quasiparticles near the Fermi level in p(i| , pl| |. However, 
there the studies were mainly addressed to the properties 
of excited states while here we will investigate only the 
properties of the ground state in the case of an attrac- 
tive interaction. We also note that the attractive case in 
one-dimension was discussed in |2^,^3| . 

To determine the characteristics of the ground state 
of the generalized Cooper problem (||) we solved numer- 
ically the Schrodinger equation. After that we rewrite 
the obtained ground state in the original lattice ba- 
sis with the help of the relation between lattice basis 
and one-particle eigenstates |n) = J2 m ^■n,m\4 > m) ■ As 
a result of this procedure we obtain the two-particle 
probability distribution _F(n.i,n.2) from which we ex- 
tract the one-particle probability /(n) = J2n 2 ^( n ii n 2) 
and the probability of interparticle distance fd(r) — 
En 2 -F( r + n 2> n 2) with r = ni - n 2 . 

Typical examples of such probability distributions are 
presented in Fig. 1. Without interaction, at given dis- 
order strength W = 2V and W = 5V both particles are 
delocalized on the lattice of given size L = 40. In pres- 
ence of interaction U = —2V the ground state remains 
delocalized for W — 2V and the probability distribution 
is rather similar to the case of U = (compare Figs. la 
and lc). On the contrary for W — 5V interaction com- 
pletely changes the ground state properties leading to a 
clear localization of both particles near each other (com- 
pare Figs. lb and le). The wavefunction is localized in 
a rather compact way and the finite size of the lattice 
definitely does not affect this localization. Fig. If shows 
that in this localized state the particles remain corre- 
lated being close to each other. This bi-particle localized 
ground state is obtained by exact diagonalization of (j^) 
where all non-diagonal interaction induced matrix ele- 
ments are taken into account. It is interesting to compare 
this solution with the mean field approximation (Cooper 
ansatz) in which only diagonal terms are taken into ac- 
count. Within the Cooper ansatz the particles occupy 
the same non-interacting orbitals and only matrix ele- 
ments Q mi , m2im ^ m ' with mi = m 2 and m 1 = m 2 are 
kept in (||) . The ground state obtained from the Cooper 
ansatz is shown in Figs.lg and lh and is clearly delo- 
calized contrarily to the strongly localized ground state 
obtained from exact diagonalization of (Q) and shown in 
Figs.le and If. In fact the ground state from the Cooper 
ansatz is more close to delocalized non-interacting eigen- 
state in Fig. lb than to the real eigenstate in Fig.le in the 
presence of interaction. The results of Fig.l definitely 
show that the attractive interaction leads to localization 
of pairs in the regime when non-interacting states are de- 



localized. This localization is not captured by the Cooper 
ansatz which neglects non-diagonal matrix elements and 
due to that misses the essential physical effect. 
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FIG. 2. Dependence of the inverse participation ratio £ of 
the TIP ground state on the disorder W/V. The full lines with 
full symbols correspond to an interaction strength U = —4V, 
the dashed ones with open symbols to U = —2V and the 
dot-dashed ones to U = 0V (+, X,*). Different symbols cor- 
respond to different linear size of the lattice L = 20 (°,+), 
L = 30 (□, x), and L = 40 (A,*). 

In order to study the ground state properties of our 
model in a more quantitative way it is convenient to com- 
pute the inverse participation ratio (IPR) £ defined as 
= (J2 a f 2 ( n )) where the brackets mark the averag- 
ing over Njj disorder realizations (typically Njj = 100). 
Physically, £ gives the number of lattice sites occupied by 
one particle in the TIP ground state. The dependence of 
the IPR £ on the disorder strength W is shown in Fig. 2 
for different strength of interaction U and different sys- 
tem sizes L. In the absence of interaction, for finite sys- 
tem sizes used in our numerical simulations (L < 40), 
the ground state is delocalized for disorder W < hV 
and it becomes localized for W > bV. On the contrary 
in the presence of interaction the TIP ground state be- 
comes localized for W > W c w 2V at U = —4V and for 
W > W c w W at U = -2V. Indeed for W < W c the IPR 
£ starts to grow significantly with the increase of the sys- 
tem size L that corresponds to pair derealization. The 
decrease of the W c value induced by the attractive inter- 
action shows that the attraction leads to localization of 
pairs inside the non-interacting delocalized phase. This 
effect is absent in the mean field approximation where the 
ground state remains well delocalized (compare Figs, le 
and lg). This phenomenon is similar to the situation in 
three-dimensional Anderson model for which the localiza- 
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tion of pairs was discussed in |T^] . As in we attribute 
this phenomenon to the increase of the effective mass 
m e ff of the pair (V oc l/m e //) that leads to a decrease 
of the critical disorder strength (W c ~ V oc l/m e //). For 
strong attraction the mass is approximately doubled so 
that the value of W c is decreased by a factor of two com- 
paring to the non-interacting case. The numerical data 
in two and three dimensions presented here and in |T^] 
are in satisfactory agreement with this estimate. 




L 

FIG. 3. IPR £ as a function of the linear size L. Dashed 
lines and open symbols are for W/V = 2 and full lines and 
full symbols are for W/V = 5V, with U/V = (o), -1 (□), 
— 2 (diamond), —4 (A). The average is done over 100 dis- 
order realizations. Here, the filling factor is v = 1/4 and 
a = L 2 /M = 2. 

The dependence of IPR £ on the system size L for 
W = 2V and W = 5V is shown in Fig. 3 for different 
values of interaction U. For L < 12 the non- interacting 
states are delocalizcd. The introduction of interaction 
decreases significantly the IPR value and for W = 5V 
the TIP ground state is localized for \U\ > 2V. On the 
contrary for W = 2V the IPR still grows with L for 
U = —2V. This behavior is qualitatively similar to the 
numerical data obtained in ps[ by projected quantum 
Monte Carlo method (see Fig. 3 there). According to JlJ] 
the pairs at quarter filling become localized at U/V k —4 
for W = hV and remain delocalizcd for W = 2V (only 
sizes L < 12 where accessible by this method). While the 
qualitative behavior is similar (compare Fig. 3 here with 
Fig. 3 in |l8|]) the quantitative difference between the two 
sets of data is definitely present. For example in our Fig. 3 
at W = 5V the states become localized approximately at 
U/V w -1.5 and not at U/V « -4 as in @. We at- 
tribute this quantitative differences to the fact that in 
fl8f up to 74 real spin fermions were present and were 



treated exactly (up to statistical errors) by the quantum 
Monte Carlo method. The presence of other fermions 
can renormalize effective strength of interaction between 
two particles. Also it can change the effective strength 
of disorder for fermions near the Fermi level. The com- 
parison between the two figures shows that the TIP ap- 
proach captures the qualitative physical properties of the 
system but quantitatively it gives different values. In a 
sense this situation is similar to the comparison between 
the Cooper approximation and the BCS theory. 
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FIG. 4. Dependence of the coupling energy A on the disor- 
der strength W for U = —2V. The linear lattice size is L — 20 
(o) and L = 40 (A); data are obtained by exact diagonaliza- 
tion of TIP model (solid lines, full symbols) and by mean field 
approximation from the Cooper ansatz (dashed lines, empty 
symbols) . 

The difference between exact diagonalization of the 
TIP Hamiltonian (0) and the mean field solution given 
by the Cooper ansatz is also clearly seen in the coupling 
energy of pair A = E g {U = 0) — E g (U). Here E g is the 
TIP ground state energy in presence of interaction U . 
For U = we have E g (U = 0) = 2E F . The dependence 
of A on the disorder strength W is shown in Fig. 4 for 
U = -2V. In the BLS phase at W > W c w W the 
coupling energy A obtained from exact diagonalization 
becomes significantly larger than the value of A given by 
the mean field approximation based on delocalized states. 
This shows that energetically the BLS phase is more fa- 
vorable than the mean field Bogolubov-dc Gcnncs solu- 
tion p3] . The physical reason for the increase of A com- 
paring to the mean field value is related to localization: 
the pairs are localized and particles remains closer to each 
other that effectively increases the coupling strength be- 
tween them |^5| . On the contrary for W < W c when 
the pairs are delocalized the exact solution gives the val- 
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ues of A which are close to the mean field value. This 
is in agreement with the Anderson theorem according to 
which mean held remains valid in the regime with weak 
disorder. 

The fact that for TIP the BLS phase is energetically 
more favorable than the mean held solution indicates that 
also at finite particle density the BLS phase will be more 
favorable. This indication is in agreement with the quan- 
tum Monte Carlo computations presented in p8| . It is 
also possible to give another argument in the favor of BLS 
phase. Suppose that the filling factor v is close to the 
critical value v c at the mobility edge of non-interacting 
particles (but v > u c ). Then it is natural that all par- 
ticles below the mobility edge are localized. Then the 
density of interacting pairs above v c is proportional to 
\v — v c \ and is relatively low for \v — v c \ <C v. In this 
regime the pairs above v c are well separated and the TIP 
approximation we discuss in this paper should be rather 
reasonable. Of course on the next step the residual inter- 
action between pairs should be taken into account [p6| . 
In this picture it is clear that the BLS phase is energet- 
ically more favorable comparing to the delocalized mean 
field solution. 

III. GROUND STATE PROPERTIES WITH 
MAGNETIC FIELD 

It is interesting to understand how the TIP properties 
in the BLS phase are affected by a magnetic field B per- 
pendicular to the two-dimensional lattice. In this case 
the one-particle Hamiltonian takes the form 

H 1= J2E n \n)(n\ + Vj^ F(n) + T») (3) 

n n 

where T(n) and T*(n) are the translation operators from 
site n to its nearest neighbours 

T(n) = T Mx (n)|n)(n + e x |+T My (n)|n)(n + e y |. (4) 

Here e x and e y are the unitary vectors on the two di- 
mensional lattice and 

^ (y) (n)=exp(-|/ r ^ (n) A.dn^ (5) 

are the magnetic translation operators along paths 
L x (n) = (n — » n + e x ) and r y (n) = (n — > n + e y ). For 
convenience we choose the Landau's gauge for the mag- 
netic field A = —n y Be x . The magnetic translation op- 
erators are then determined as Tm x = exp{2 r Ki r yn y ) and 
Tfyfy = 1 with 7 = Due to the periodic boundary 
conditions the effective topology of the two-dimensional 
lattice is that of a torus with a transversal and a lon- 
gitudinal radius Rt = Ri = L. This topology implies 
flux quantization on the lattice so that 7 = jj with 



m G [0, L — 1]. Then the one-particle Hamiltonian (||) 
can be explicitly written as 

ff 1 =^ J B n |n)(n| 

+ ^(e 2 ^|n)(n + e x | + |n)(n + e y |) 

n 

+ E (V 2 ™ 7 ™" |n + e x ) <n| + |n + e y } <n|) . (6) 

n 




FIG. 5. Ground state probability distributions of two inter- 
acting particles in a lattice of linear size L — 40, with disorder 
strength W = 3V , and interaction strength U — —2V; the left 
column shows the probability distribution /(n) and the right 
one the interparticle probability distribution /d(r). The dif- 
ferent rows presents the TIP ground state for different value 
of the magnetic flux 7: top row 7 = 0, middle row 7 = ^j, 
and bottom row 7 = 

We study now the ground state of the TIP Hamiltonian 
constructed with this new one-particle Hamiltonian H\ 
(H). This Hamiltonian is written in the preferential ba- 
sis of non-interacting eigenstates of (||) that leads to the 
Schrodinger equation of the form (Q). As in the case with 
B = we use the probability distributions /(n), /d(r) 
and the IPR £ to depict the ground state properties of the 
TIP problem in presence of a magnetic field. Thus Fig.|5| 
represents the TIP probability distributions for a system 
of linear size L — 40 and fixed disorder and interaction 
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strengths (W = 3V, U = — 2V). The data are shown for 
different magnetic flux ratios 7 = 0, 7 = ^ and 7 = Jj. 
At 7 = the one particle probability is well localized 
by interaction (first row of Fig.^J). However, with the 
increase of magnetic flux 7 the localization is destroyed. 
The data of Fig|| suggest that there exists a critical mag- 
netic flux 7 C below which the TIP pairs remain localized 
(7 < 7c, middle row of Fig.^J) and above which pairs be- 
come totally delocalized (7 > j c , bottom row of Fig.||). 
At the same time for 7 > "f c the interparticle distance 
probability distribution /<j(r) is less peaked. Hence for 
7 > 7 C the size of the pair is significantly increased. 
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FIG. 6. Dependence of the inverse participating ratio £ of 
the TIP ground state on the disorder W/V and magnetic flux 
7. The interaction strength for the main figure is U = — 2V . 
Empty/full symbols correspond to 7 = 0/7 = 0.1, with differ- 
ent linear lattice size L = 20 (o) and L — 40 (A). Inset shows 
the phase diagram in the plane of W c /V and 7; the full line 
corresponds to U = —2V and the dashed line to U = —4V; 
BLS phase is at W > W c . 

The ground state properties can be studied in a more 
quantitative way with the help of the 1PR £ defined 
above. FigJ| represents the dependence of £ on the dis- 
order strength W for a fixed interaction U — —2V. Data 
are shown for different values of magnetic flux 7. They 
clearly show that the introduction of magnetic field leads 
to an increase of £ at a fixed value of W. Thus the 
magnetic field enhances the derealization of particles 
for 7 > 7c (U) and W < W C (U). On the contrary for 
7 < 7c(c0 and W > W C (U) the variation of £ with lattice 
size L is weak and usually here there is a small decrease 
of £ with increase of L (see Fig.^). The derealization 
transition can be determined as the point where £ is in- 
dependent of the lattice size (crossing point). An approx- 
imate phase diagram in the plane (W, 7) obtained in this 



way is shown in the inset of Fig.[| With the increase of 
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FIG. 7. Dependence of the inverse participation ratio £ of 
the TIP ground state on the magnetic flux 7 for two lattice 
sizes L — 20 (•) and L — 40 (triangle). Here the interaction 
strength is U = ~2V and the disorder strength is W = 3V. 
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FIG. 8. Pair coupling energy A as a function of magnetic 
flux 7 for L = 20 (•) and L = 40 (triangle). The interac- 
tion and disorder strengths are as in Fig.|7j: U = —2V and 
W = 3V. 

interaction the localized phase penetrates deeper in the 
region of weak disorder. The fact that a magnetic field 
can delocalize pairs inside the BLS phase is also illus- 
trated in Fig.fT] where the IPR £ is enormously increased 
by the magnetic flux. This result is in agreement with a 
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general fact known for non-interacting particles that the 
localization length is increased by a magnetic field Jl9| . 
While the delocalization induced by a magnetic field is 
clearly illustrated by the FigsJ| and [?] it is rather difficult 
to determine numerically the properties of pairs in the 
delocalized phase. In this phase an effective pair size be- 
comes too large to investigate it numerically. The ques- 
tion if the superconductivity survives or if a magnetic 
field drives the system to a metallic regime is difficult to 
answer in the frame of our numerical approach. 

For a better understanding of both localized and de- 
localized phases we studied the dependence of the pair 
coupling energy A on the strength of magnetic field. In 
the standard Cooper problem A is related to the BCS 
gap and determines the Cooper pair size l pa ir oc I/A. 
Fig.|§| shows the dependence of A on the magnetic field 
for different lattice sizes and for the case U = — 2V and 
W = 3V, already used in the Fig.fj]. In the localized 
regime with 7 < 7 C 0.04 the value of A varies weakly 
with the growth of L. In contrast, for 7 > j c its value 
decreases in 2 - 3 times. This can be considered as an 
indication that the superconductivity is significantly sup- 
pressed by magnetic field for 7 > j c . However, a signifi- 
cant increase of L is required to investigate the properties 
of pairs in the delocalized regime. 



IV. CONCLUSION 

The present studies show that in the presence of disor- 
der the attractive Hubbard interaction leads to localiza- 
tion of pairs and appearance of phase with bi-particle lo- 
calized states which is located inside the non-interacting 
metallic regime. This BLS phase cannot be obtained in 
the mean-field approximation. It is shown that it can be 
destroyed by the introduction of a magnetic field which 
drives the system to delocalization. In the BLS phase 
the pair coupling energy is much larger than the value 
obtained in the mean-field approximation. This indicates 
that the BLS phase is energetically more preferable com- 
paring to the mean-field solution. The results obtained 
for two particles (one pair) are in qualitative agreement 
with the recent results obtained with the quantum Monte 
Carlo in f|. 
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